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ABSTRACT 

In this article we extraet soft distribution functions for Drell-Yan and Higgs production processes 
using mass faetorisation theorem and the perturbative results that are known upto three loop level. 
We find that they are maximally non-abelien. We show that these functions satisfy Sudakov type 
integro differential equations. The formal solutions to sueh equations and also to the mass faetori¬ 
sation kernel upto four loop level are presented. Using the soft distribution function extracted from 
Drell-Yan production, we show how the soft plus virtual eross seetion for the Higgs produetion 
ean be obtained. We determine the threshold resummation exponents upto three loop using the 
soft distribution funetion. 



The Drell-Yan(DY) production of di-leptons and Higgs boson production play crucial role in 
the hadronic colliders. The di-lepton production can not only serve as a luminosity monitor but 
also provide vital information on physics beyond standard model at present collider Tevatron at 
Fermi-Lab and future Large Hadron Collider (LHC) which is going to be up at CERN in few 
years. Higgs production at such colliders will establish the Standard Model(SM) as well as beyond 
SM Higgs [1,2]. From the theoretical side, the DY production of di-leptons and Higgs boson 
production are known upto Next to Next to leading order(NNLO) level in QCD. For DY at NFO 
level, see [3] and for the Higgs production at NFO level, see [4-6]. The NNFO contribution to 
DY can be found in [7-9]. Beyond NFO, the Higgs production cross sections are known only 
in the large top quark mass limit. For the NNFO soft plus virtual part of the Higgs production, 
see [10,11] and the full NNFO for the Higgs production can be found in [12-14]. Apart from 
these fixed order results, the resummation programs for the threshold corrections to both DY and 
Higgs productions have also been very successful [15,16]. For next to next to leading logarithmic 
(NNFF) resummation, see [17,18]. Due to several important results at three loop level that are 
available in recent times [19]- [24], the resummation upto N^LL has also become reality [25-27]. 

With all these new results in both fixed order as well as resummed calculations, one is now 
able to unravel the interesting structures in the perturbative results (for example: [34-36]). Along 
this line, in this paper, we extract the soft distribution functions of Drell-Yan and Higgs produc¬ 
tion cross sections in perturbative QCD and show that they do not depend on the process under 
consideration. By that we mean that the soft distribution function of Drell-Yan production can be 
got entirely from the Higgs production by a simple multiplication of the colour factor Cp jC a- We 
prove this for the pole parts upto three loop level and for the finite part we could show only to 
those terms that are not proportional to 5(1 —z) because the three loop finite part proportional to 
5(1 — z) is not available yet and can be obtained only from the explicit fixed order computation 
of bremsstrahlung contribution. The extraction of the soft distribution function is achieved with 
the help of mass factorisation theorem supplemented by the recent developments in the compu¬ 
tation of three loop anomalous dimensions, three loop form factors of quark and gluon operators 
and two loop bremsstrahlung contributions to Drell-Yan and Higgs productions. We discuss the 
consequences of our observation in the determination of soft plus virtual cross sections and the 
threshold resummation exponents. A brief account on the soft and jet distribution functions and 
the resummation exponents relevant for deep inelastic scattering (DIS) is given. 


We start by writing the partonic cross section as 




A) 





as,Q^,l?) p5(l-z)(g)Le 





l = q:g ( 1 ) 

with the normalisation, = 5( 1 — z). The symbol sv means that we restrict to only the soft and 
virtual contributions to the partonic cross sections 6^''. In the above equation we have introduced 
a "L ordered exponential" which has the following expansion: 
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The function /(z) is a distribution of the kind 5( 1 — z) and “Di, where 


‘Di = 


W(l-z) 

. (1-z) . 


/ = 0,1,- 


+ 


(3) 


and the symbol ® means the Mellin convolution. The letters q and g stand for Drell-Yan(DY) and 
Higgs(H) productions respectively. q^{= —Q^) is the invariant mass of the final state (di-lepton pair 
in the case of DY and single Higgs boson for the Higgs production), z is the scaling variable defined 
as the ratio of q^ over s, where s is the center of mass of the partonic system. F\ds, are the 

form factors that enter in the Drell-Yan(for I = q) and Higgs(for I = g) production cross sections. 
The functions {ds^q^^z) are called the soft distribution functions. The unrenormalised(bare) 
strong coupling constant ds is defined as 


av = 


'^2 

si 

1671^ 


(4) 


where gs is the strong coupling constant which is dimensionless in n = 4 4- £, with n being the 
number of space time dimensions. The scale /r comes from the dimensional regularisation in order 
to make the bare coupling constant gs dimensionless in n dimensions. 

The bare coupling constant ds is related to renormalised one by the following relation: 

Stds = Z{n\)as{i?R) (5) 


The scale fUR is the renormalisation scale at which the renormalised strong coupling constant 
is defined. 


Se = exp 



[Y£-ln47t]| 


( 6 ) 


is the spherical factor characteristic of n-dimensional regularisation. 

The fact that ds is independent of the choice of pR leads to the following renormalisation group 
equation (RGB) for the coupling constant: 


Fr- 


d I n as{pl 

dpi 


= - + 


2 asipl) 


^(asipl)) 


where 


= -asipl) pI 


dlnZ{pl) 


dpi 


i=0 


The solution to the above equation is given by 


ApI) = 1 + asipl) ^ ipl) 




3£ 


(7) 


( 8 ) 


(9) 
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The renormalisation eonstant relates the bare eoupling eonstant as to the renormalised one 
as{ia\) through the eqn.©. The eoeffieients (3o and Pi are 

n 11 4 

Po = ~ 3^^”/ 

Pi = — — —TfHfCA (10) 

where the eolor faetors for SU (N) QCD are given by 

N^-l 1 

Ca=N, Cf = -^, Tf = - (11) 

and Hf h the number of active flavours. In the case of the Higgs production, the number of active 
flavours is five because the top degrees of freedom is integrated out in the large rtitop limit. 

The factors the overall operator renormalisation constants. For the vector 

current Z‘i{ds,iJp,iJ^) = 1, but the gluon operator gets overall renormalisation [28] given by 


Z^ds,/jl,p^,e) 
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, -'3 
+a^ 



2('-2PoP,)+?^ 


( 12 ) 


The bare form factors F\ds, (before performing overall renormalisation) of both fermionic 

and gluonic operators satisfy the following integro differential equation that follows from the gauge 
as well as renormalisation group invariances [29-32]. In dimensional regularisation, 


d 








+ G^ 




(13) 


where contains all the poles in £. On the other hand, collects rest of the terms that are finite 
as £ becomes zero. In other words contains only non-negative powers of £. Since is RG 
invariant, we find 
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The quantities are the standard cusp anomalous dimensions and they are expanded in powers of 
renormalised strong coupling constant asifip) as 




a' 


i=l 


(15) 
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The total derivative is given by 


^^R 


T 2 2 
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The RGB of can be solved in powers of bare eoupling eonstant ds as 




.P 

,2 \ 


Sh 


where, 

K‘d(e) = 

K‘^^{e) = ^|^2(3oAi^+-|^-A2^ 

/^^A(g) = _ -p _ ^ _ -PoP^Aj — 6P^2^ 

P2A1 + P1A2 + SPoAj^ e ( ~ 2"^^) 
Similarly RGB for ean also be solved and the solution is found to be 


= G‘ [asiO^), 1 ,£) + {usiX^nl)) 

The integral in the above equation ean be performed and it is found to be 
^^A\as(k^nl)) = 

The finite funetion G\as{Q^), 1,£) can also be expanded in powers of as{Q^) as 




G'(a,(e2),i,e) = £4(e2)Gf(£) 

i=\ 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


(21) 
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After substituting these solutions in the eqn. (fT3t and performing the final integration, we obtain 
the following solution 

oo 

lnF'(a„e2,;/2 £) = ^4 

i=i 



^ Sc 


£ -^F 


/,(0 


(e) 


( 22 ) 


where 

^ ^ PoG5(e)^ + —G2(e) 

= i(-|p5A!)+^0|5,A!H-|poA'2 + l|55Gi(E)) 

4(-|A5->G!(e)-ip„G'2(e))-Ei(iG5(e)) 

+i3 (^■^P 2 ^i + ^Pi 2^2 + ;^PoA 3 + -PoPiG{(e) + 3PoG2(e)^ 

The above result is in agreement with [21], whieh was evaluated using various algorithms designed 
for solving nested sums. The cusp anomalous dimensions and G{(£) are known upto order 
a^. The cusp anomalous dimensions are maximally non-abelien and hence satisfy the following 
relation: 





(24) 


The coefficients G((£) can be found for both I = q and 7 = g in [22] to the required accuracy 
in £. They satisfy 

oo 

G\ (£) = 1 {b{ — 5/,g|3o) + /i + 52 

k=\ 

OO 

a‘M = 2(B'2-28,,jP,)3-/l-2P„*,'’‘ + £eV'‘ 

k=\ 


5 


(25) 


G'^{e) = 2(55-38,,j|32)+/j-2|5i«,'-'-2|3o(«2'’‘+2|3o«,'’") 


+ 




k=l 


The constants are also known upto order thanks to the recent computation of three loop 
anomalous dimensions/splitting functions [19,20]. 

The constants fj are analogous to the cusp anomalous dimensions that enter the form factors. 
It was first noticed in [33] that the single pole (in £) of the logarithm of form factors upto two loop 
level (a^) can be predicted due the presence of these constants fj because they are found to be 
maximally non-abelien obeying the relation 

ff = (26) 

similar to Aj. In [22], this relation has been found to hold even at the three loop level. 

The partonic cross sections is UV finite after the coupling constant and over¬ 

all operator renormalisations are performed using Z(ju|) and Z^(/r|). But they still require mass 
factorisation in order to remove the collinear divergences: 

with /jp being the factorisation scale. The resulting coefficient functions are finite 

and free of collinear singularities. 




/=l 




(28) 


The coefficient functions for / = 1,2 are known(see [4] to [14]). The partial result for 
(i.e., all £>,■ except 5(1 —z) are known for i = 3) is also available(see [25]). 

The kernel r(z,/r|., e) satisfies the following renormalisation group equation: 

^F^r(z,A/^,£) = {z,i4) Z>r (z,A'F.e) (29) 

The P{ztIJ^f) known Altarelli-Parisi splitting functions(matrix valued) known upto three 

loop level [19,20]: 

CX3 

P(z,A/^) = £4(A/^)p('"i^(z) (30) 

i=\ 


The diagonal terms of splitting functions P^^^ (z) have the following structure 


pri\z)=2 


^/+lS(l -z) +A\^-^Pq 


+ P, 


reg,II 


(z) 


(31) 
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where are regular when the argument takes the kinematic limit(here z —1). The RGB of the 
kernel can be solved in dimensional regularisation in powers of strong coupling constant. Since we 
are interested only in the soft plus virtual part of the cross section, only the diagonal parts of the 
kernels contribute. In the MS scheme, the kernel contains only poles in £. Expanding the kernel in 
powers of bare coupling ds. 


r(z,A'F,e) =§(l-z) + L 


i=l 


2 \ 'T 


we can solve the RGE for the kernel. The solutions in the MS scheme are given by 


dfU.e) 




(32) 




1 /I 


e\2 


r!,”(2,E) 


r!f(z.e) 


4f¥Po4“’w-Po<’(z)®d"’w 




(z) ® P^\z) ® p,f- (z)^+^ (if,',"’ (z) 0 f,'/’ (z) 


-i|5i/i,"’(z) - i|5o/i,‘’(z)) + i(if,*i”(z) 


■i|3„4"’(z) 0P,<,">(z) ®P<,"’(z) + E|35 /^,»'(z) ®f<,"’(z) 




(if,',"’ (z) ® f,',"’ (z) ® f,','’(z) - ipif,',"’ (z) ® f,',"’(z) 


■Tp„p|,») (j) + (j) + 3|32 /^;>(z) 


+ - 


I (if,',"’ (z) ® fj? (z) + if,';’ (z)»pI!>(z)- ipz/i,"’ (z) 


■iPi4”W-iPo4’”w)+i(if,';’W 


(33) 



It is now straightforward to obtain the soft distribution functions z) from the available 

results known up to three loop level for the form factors the kernels Tjj and the coefficient 
functions A;^’'(the 5(1 —z) function part of is still unknown). The fact that A)'' are finite in 
the limit £ —;> 0 implies that the soft distribution functions have pole structure in £ similar to that of 
and F//. Also, satisfy the renormalisation group equation: 


=0 (34) 

dn% 

From the above observations, it is natural to expect that the soft distribution functions also sat¬ 
isfy Sudakov type integro differential equation that the form factors F^{Q^) satisfy(see eqn. lfTSt i. 
Hence, 




_d_ 

dq^ 






ds,^,Z,£^ + G^ %^,^,z,£ 


(35) 


where again K ^ contains all the singular terms and G ^ are finite functions of £. The renormalisa¬ 
tion group invariance leads to 


fI 


_d_ 

d^l 


K 





-A\as{ijl))d{l-z) 


Fr^G^ (ds,^,^,z,£) =A\asiijl))d{l-z) (36) 

dFR \ Fr F^ J 

If ^>^(dj,<?^,^^,z,£) have to contain the right poles to cancel the poles coming from and F// 
in order to make A)^ finite, then A ^ have to satisfy 

a' ^-A^ (37) 


The above relation along with the renormalisation group invariance implies that 

( d ^,^,^, z ,£) = (^s{f\),^,z,£ 


= g' {as{q^), l,z,£) -5(1 -z) ^^2 


(38) 


Now it is now straight forward to determine all G \as{q^)^ l,z,£) from the available informa¬ 
tions. The functions G {as{q^), l,z,£) can be expanding in powers of as{q'^) as 

CX3 

G\as{q^), l,z,£) = Y.<{q^) g!{ z,£) (39) 

i=\ 
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The solution to the eqn (l35b ean be obtained in the way we obtained lnF^(2^). Expanding the soft 
distribution funetions in powers of bare eoupling ds as 


£ 

oo / 2 \ ^ 2 

d>' {d,,q\^i\z,e) = (% ) 5^ 6'’«(z,£) 

,=1 VA' / 

we find the solution: 

4>'’«(z,£) = 2;’('^(£) (a^ ^ -5(1 -z) A^, G\e) ^G\z,e) 


(40) 


(41) 


The finite funetions G, (z, £) ean be obtained using the mass faetorisation formula by demanding 
the finiteness of the eoeffieient funetions . The RG invarianee of theses soft funetions and the 
simple rescaling q ^ {I —z)q imply that the following expansion is also the solution to the integro 
differential equation: 


^‘{d,,q^,p^,z,e) = 4)'(d„^^(l-z)^/,£ 


2 .2 


= 

i=\ 




2X'i 




I £ 

1 -z 




(42) 


where 


^fi(0(e) = Lp^'\t) ( a‘ -A‘,G\z) Q^z) 


(43) 


The z independent constants g (£) in $ ^40 (^e) can be obtained using the form factors, mass factori¬ 
sation kernels and coefficient functions expanded in powers of £ to the desired accuracy. 

This is achieved by comparing the poles as well as non-pole terms in £ of $ ^40 (^e) with those com¬ 
ing from the form factors, overall renormalisation constants and splitting functions and the lower 
order We find 


Q \ (£) 


kTZ I,(k) 




k=\ 




Zvfi(l) 


a— 
'2 


k=l 


= -/3-2Pi£;,'’‘'’-2|3o((;2'’"’+2|3o(;,'’'’’) + £e‘£;3'’'‘’ 

^ ' k=l 


(44) 


with 


—7,(1) — (1) 

fi 1 = <-/ fi 1 
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469 


^2 + 4^1 



+Cjnf 



70^ 32^ 

+ -^^2+ -^^3 


where Q = Cf for I = q(DY) and C/ = Ca for I = ^(Higgs). 

Using such compensating $^40 (^e) and the following expansion, 




■D, 


=0 J- 


(45) 


(46) 


we obtain {z,t) upto three loop level. We find that the finite functions g/(z,£) have the fol¬ 
lowing decomposition in terms of cusp anomalous dimension A\ and // that appear in the form 
factors: 

oo 

Gi = —/[ 5(1 — z)-|-2Aj ©0+^ e^g/’*(z) 

k=l 

oo 

G 2 = —/2 5(1 — z)-t-2A2 — 2Po ,?/’^(z) + 52 

k=\ 

G 3 = -/g 5(1 -z)+ 2A^ ®0-2pi g/’^(z)-2po(^2^’^(z)+2|3og/’^(z)) 

00 

+ (47) 

yt=l 

where 

= q(^8®i- 3^25(1-z)^ 

= C/^-3^2®0 4-42)2+2^38(1-z)^ 
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- 1,3 
Si = 


Cl( -3^2®1 + -z) 
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1616 242^ \ /1072 
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27 
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/2428 469^ 176^ A ^ //224 44^ \ 
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2 1210 
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^ +Cinf 


656 140^ 220^ , 
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3 ^2^2 


/7744\ , ^ ^ 

+ ( —^ j®3 1 +CiCAnf 


//96482 2452 1264 2 6536 
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9 ^ 3)^0 


72008 9056^ 448^ \ / 22400 320^ \ / 2816\ 


+C/n^ 


4480 1520^ 416^ \ /4192 736 

+ ^^^2 + ^^3]‘Do + 
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1- 


81 


^2 ) 


+ - 


1600\ 


27 


256 \ 
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96.0 304, 


J®2+ +C/C/rn/|^|^—^-60^2-^^2-^^3j®0 


+ -220+ 192^3 j ®1 + (^64j ® 2 j +5gf95(l -z) 


(48) 


In the above equation 5,^2 ,5,^3 are not known because the full fixed order N^LO computation 
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for the soft part of the eross seetion is not available yet. 

With the available informations (ignoring we find that the soft distribution func¬ 

tions for DY and Higgs productions are maximally non-abelien: 

^ (49) 

upto three loop level. At the cross section level(A)^), this property does not show up because of 
the form factors which do not have this property. The overall factors Cp and Ca ordinate from the 
leading order contributions to the soft distribution functions. Hence if you factor out this colour 
factor (Cf for the DY and Ca for the Higgs) we find that the soft distribution functions are universal. 

=Cj^{ds,q^,f/,Z:e) (50) 


The universality of the soft distribution functions can be understood if you notice that the soft part 
of the cross section is always independent of the spin, colour, flavour or any other quantum numbers 
after factoring out the bom level cross section. It depends only on the gauge interaction, here it is 
SU (N). This universal property can be utilised to compute soft part of the any new cross section 
where incoming particles carry any spin,colour,flavour or other quantum numbers. For example, 
if we know $(dj,^^,/r^,z,£) extracted from the Drell-Yan production results, we can predict the 
A^^(z,< 5 r^,/r|,/r|.) for the Higgs production using mass factorisation formula provided we know the 
gluon form factor F^{ds, and the overall renormalisation constant £). 

The soft plus virtual part of the cross section (A)^(z,( 2 '^,/r|,/r^)) using mass factorisation for¬ 
mula is found to be 


^7 {z,q^ 



(51) 


where is a finite distribution. The dependence comes from the coupling 

constant and operator renormalisation: 


^\z,q\A,A,^) = (\n(z\asA,F\^)] +\n\F\as,Q^^i^^)nh{\-z) 


F2Ci^{ds,q^,1?,z,t) - 2 C\nTii{ds,n^,i?p,z,t) (52) 

In the above equation "Cln" means the "convolution ordered" logarithm. All the products of dis¬ 
tributions in the logarithmic expansion are understood as Mellin convolutions. The distribution 
'F^(z,^^,/r|,/r|.,£) is regular as £ -^ 0. The soft plus virtual cross section can be obtained by 
expanding '^^(z, £ = 0) as 

(53) 

i=\ 

where we have set ^ir = np and expressing as{i^p) in terms of asi/J^) is straightforward. We find 
that the cross sections l^^f{z,q^,ij7p) can be obtained using 

^7'^^\z,q^,^7p) = C/5(l-z) 
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where 

VJ / /,( 1 ) 


VJ / ^ i ( 2 ) 


VJ / ^,( 3 ) 










C'” (z,rA) = f (z, + f '■“* (z, 0 f (z, ) 


+if'’">(z,?^#/|) 0f'’<‘>(z,9^^/?■) 0f'■“'(z,«^^/?•) 


(54) 


-2Po8;,s + 2B')5(1 


+ 2 C ;(7 ^^5(1 —z) + ^4A^ 


In + (^3^24^1 +2g/’ 
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6^2^0-81 ~6^2Po^/,g + 2 Po,?i ’ +,?2 + 3 ^ 24^2 + 2 ( 3 oC/^/ ^+C/^2^ ^ 


+ ( -4(3i5/,^-2porf425i-3C2Mi-2|3oC/^;'Mln 


7t(i) 





+ ( — P 081 + Po5/,g) In 


,2 / ^ 



5(l-z) + 


^ ( 1 ) o 2 -/ I o 4 / 1 *. ^ 


-4|3oC/^i^ ^-2/^+2A^2ln 



-poA^ln^ ( 



®o + 
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4A'2-4Miln 
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-4Mi 


®2 


30C2PoPi8,,s + nCzPoBi - nfePfe/'"’ + 6C2 PiB'i + ipo«2+ jPsg,'’*’’ 


A-^Pl.?! ^ + 2*^3 ^ +3^2A3 — 3^2P^1 + 2 PoC/^2 ^ 1 


(/ or2o2x/ 


I 8o2y-._^ (3) 


+-PiC/^/ ^ + -C/(^3^ ^+ 6^2 Po/ 2 + ^12^2PoS/,g — 6(325/,g+ 253 — 12^2Po8i 


-2(.„gl-‘-''’ - 4Pfe/’'2> - 2Pi«- e&PoAj - SfcPiA'i - 

2 



\ ) + ( 5popi8,,, - 2PoB^ - P.B'i + 2p52;/'“' 
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+ (|P6B'l-|pfcV'>g’''^ 
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-4poC/^2^^^ - - 4|3iC/(^ - 2/1 + (^4po/l + S^lQg 

+2^0 f “ 2(3 oA 2 - PiAj^ In^ f + f -P^i^ In^ (\ 


-( 1 ) 
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4) ' V "" V WrJ ' V3- 
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8 P 0/2 + 16(3qC/^ ^ + 4 A 3 + ( — §^ 2^2 ~ 4PiAj) In ( 

Pf 
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.(4P^!)ln^(| 


©1 + 


8M2-4Mi + (8p^i)ln(\ 

-.f^F 
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®2 


+ 




®3 


(55) 


Notice that and are eompletely known. Using our and we eould 

sueeessfully reproduee soft plus virtual eross seetion (see [ 10 , 11 ]) (i = 1 , 2 ) of the 

Higgs produetion from that of DY [ 8 ] and viee versa. To eompute 
we need to know and ( 7 /’*'^^ 


{z,q^,IJf) (equivalently 
r 3 either from DY or Higgs produetion beeause these 
eonstants are maximally non-abelien. Notiee that these eonstants appear only in the eoeffieient 
of 5(1 — z) part of 'T Sinee the eoeffieients of ©,(/ = 0,1,2,3) in 'T do not depend on 
these unknown eonstants [ 73 ^’*'^^ and g 2 ^'^\ we ean prediet these eoeffieients (say for the Higgs 
produetion) by using the universal soft distribution funetion extraeted from a proeess(say DY), the 
three loop form faetors and the renormalisation eonstants. Our predietion agrees with the partial 
N^LO soft plus virtual results [25] for DY and Higgs produetions. 

With these available informations one ean also determine all the quantities upto N^LL level 
in the threshold resummation. To do this, we first reeolleet that the soft distribution funetion 
is renormalisation group invariant. Its UV divergenee ean be removed by the eoupling eonstant 
renormalisation. This introduees a renormalisation seale i^r whieh is arbitrary to all orders in 
perturbation theory. In order to eompute various quantities in the threshold resummation formula 
from the soft distribution funetion, we ehoose = i^f- With this ehoiee, one ean express the soft 
distribution funetion as a sum of pole and finite parts in £ as £ -^ 0 , that is 

t q^ 1 

^s\I^F) 1 2 ’ 




q 


Mf 


^ pole 




+ ^fin 


q 


A'f 


(56) 


With this deeomposition, it is now straightforward to identify the finite part with the threshold 
resummation formula as 
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(57) 

where the subscript S in Hg indicates that it comes from the soft part of the cross section. The 
remaining contribution comes from the form factor. D\as{q^{l —z)^)) can be expanded in powers 
of bare coupling constant as follows: 

oo 2 2 i- 

The finiteness of after coupling constant renormalisation demands that it satisfies the following 
expansion in £: 

CXD 

£>/,W(e)= £ (59) 

Using RG invariance, the coefficients of negative powers of £ can be evaluated as 




-4po -4p2 - (3i 


We find that for non-negative powers of £(j > 0), 


dr = 

Using the above equations, we find explicitly 
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g,m 
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= ci 
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(62) 


^ f i ) 

The coefficients d - for the DY can be obtained using 


j Ca ^ 


(63) 


because the soft distributions functions are maximally non-abelien. Also, the coefficients of 
a\{q^)‘DQ in the soft distribution function are related to the coefficients D( that appear in 
threshold resummation formula. Hence it is straightforward to obtain from the soft distribution 
function We find that D\ are related to and hence g/’* as follows: 


D\ = 


4'’(i)=2g/’i 
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(64) 
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From the available informations on three loop results, we find that the following result holds 

D\ = l~g\{t = Q) z = 1,2,3 (65) 

The faet that D\ ean be expressed entirely in terms of g/(z,e) (i.e., in terms of or = 0)) 
whieh are maximally non-abelien, implies that D\ are also maximally non-abelien. Henee, 


' Ca ‘ 


( 66 ) 


with 
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The above results are in agreement with [25-27]. We also find the resummation exponents D\ can 
be extraeted by using the following relations: 


D = A 


5V,(1) 


®0 


D‘, = 


A^v,( 2)_ l^.v,(l)^^.v,(l) 


©0 


^3 


Af (3) - Af’)!) ® 1 Af ® ® Af 


©0 


( 68 ) 


AV i l') "7 "7 O 

where A^ in the above are eomputed at the seale /af- = iaf( = q . From the following eonvolution 
identity [37] upto irrelevant regular terms(denoted by ■ ■ ■) 


^+,/+l 

‘Di^'Dj =dijd{l-z)+ Cijj-Di-l - (69) 

1=0 

it is interesting to notiee that in order to obtain D^, it is suffleient to know the eoeffieients of all 2 )/ 
(( = Imax to 0 ) (that means, we need not know the information on the eoeffieient of 5( 1 — z) funetion 
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and the regular part of and the eomplete soft information of ^\i.e., the eoeffieients of 
all and 5(1 —z) are needed). 

Finally, the eoeffieient of 5(1 — z) in the resummation formula ean be obtained from by 
defining the eoupling eonstant at the seale The result is 


^F/ 1=1 


(70) 


where 


and 


^ 1,1 


^ 1,2 
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(72) 


The remaining contribution to the exponent comes from the the finite part of form factor. 

We conclude our discussion on this subject with a brief discussion on the corresponding soft 
as well as jet distribution functions that appear in deep inelastic scattering. The soft plus virtual 
coefficient function that appear in the hadronic structure function F 2 can be expressed 

as 


ClncY, 2 iQ\z) - +ln|F'(4,Ejp)5(1-z) 


+2d>5y(dj, 2^,^^,z,e) — C lnr//(dj,^^,^^,z,e) (73) 

where d>^y(dj, Q^,fj^,z,£) is sum of soft and jet distribution functions. We find that this soft plus jet 
distribution function also satisfies Sudakov type integro differential equation (see eqn. dSlD which 
can be solved in the same way we solved soft distribution functions. We find that this soft plus jet 
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distribution function can be expressed as 


i=l 


A' 


'■I 


= i4(^^] Sk 


l £ 


2 (1-z) 


I ^’W(£) 


where 


|/,( 0 (£) = ^ -A^,G\£) 


(74) 


(75) 


We find that the constants ((^(e) have the following expansion in terms of Bj, // and the £ depen¬ 
dent part of lower order coefficient functions. 

CO 
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The z independent constants are computed using the coefficient functions Q) known 
upto three loop level [23,38]. Recollect that the three loop form factors were obtained from these 
coefficient functions by demanding the finiteness of the partonic cross sections after mass factori¬ 
sation and also notice that the method used there is very different from the method presented in 
this paper. We obtain 

gf’W = Cf (I -31^2" 


:9,(2) 
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Using the following decomposition, 
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it is now straightforward to identify the finite part ^sj fin the DIS threshold resummation 
formula as 


dzz 


,/v-i 


^SJ,fin 




■1 
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'0 
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^DIS [^s[Q ( 1 ~ z ) 
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Using the above equation, we find that the resummation eonstants 5^/5; satisfy the following 
relation 


Blis4 = ^fh = 0) /= 1,2,3 (80) 

The resulting -s agree with those given in [23]. 

To summarise, we have extracted the soft distribution function using mass factorisation for¬ 
mula for both Drell-Yan as well as Higgs productions within the framework of perturbative QCD. 
This is possible now thanks to various three loop results available for the form factors and splitting 
functions. The is known completely upto two loop level. Except the soft bremsstrahlung contri¬ 
butions proportional to 5(1 —z) (at three loop level), all the other soft terms(!D/) are known for the 
soft distribution functions upto three loop level. We have also shown that the soft distribution 
functions satisfy Sudakov type integro-differential equation that the quark and gluon form factors 
satisfy. We found that they are process independent. In other words, knowing the soft distribution 
function of the Drell-Yan process, one can obtain the same for the Higgs production by simply 
multiplying the colour factor combination Ca/C/?. Hence, unlike the cross sections A^'", these soft 
distribution functions are maximally non-abelien. Notice that the form factors have this property 
only at the pole level. Using the universal soft distribution function extracted from DY production, 
we could reproduce the soft plus virtual cross section for the Higgs production using the gluon form 
factor and its renormalisation constant. We have also shown how these soft distribution functions 
are related to threshold resummation formula. We have extracted various coefficients appearing in 
threshold resummation formula upto three loop level using the soft distribution function. We have 
also discussed the DIS soft and jet distribution functions in the present context. 

Acknowledgments: The author would like to thank S. Moch and A. Vogt for providing the 
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